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Calcutation of kaon matrix elements in quenched domain-wall QCD with 
DBW2 gauge action 

J. Noaki a for RBC Collabolation* 

a RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973-5000, USA 

We give a progress report of our new a~ « 3 GeV quenched calculation of kaon matrix elements with domain- 
wall fermion and DBW2 gauge action. Our smaller lattice spacing allows us to address the effect of charmed 
quark on the lattice. We show preliminary results of Bk renormalized non-perturbatively and K — * n matrix 
elements. 



1. Introduction 

In order to treat the QCD effect in non-leptonic 
kaon decays non-perturbatively, it is important 
for the matrix elements of the local operators to 
be calculated on the lattice. A couple of years 
ago, CP-PACS and RBC Collaboration jG] cal- 
culated all of the matrix elements for the interac- 
tion of K — > 7T7T decay: Hw = J2i @i(fi)Qi- Using 
the domain-wall fermion formalism to realize the 
chiral symmetry required in this calculation, they 
reported small and negative values of e'/e in con- 
flict with the experimental result. In these works, 
however, there are several uncontrolled system- 
atic errors coming from 1) the small, but non-zero 
breaking of chiral symmetry, 2) finite lattice spac- 
ing, 3) the perturbative treatment of the charmed 
quark in the matrix elements, 4) quenching effect, 
and 5) K — > 7T7t matrix elements obtained from 
K — > 7r and K — > by using lowest order chiral 
perturbation theory 

In order to examine the first three of above 
systematic errors extensively, we are performing 
a quenched simulation with domain-wall fermion 
and the DBW2 gauge action which improves the 
chiral symmetry on the lattice 0]. The degree 
of chiral symmetry breaking is decreased by a 
factor 1/10 compared with the previous work of 
RBC Collaboration. In addition, the effect of the 
charmed quark on the lattice can be examined as 
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well as the scaling violation. 

In the rest of this article, we present the con- 
tents of the numerical simulation and report pre- 
liminary results of kaon B-parameter Bk and 
the matrix elements which numerically dominate 
AI =1/2 rule and e'/e. Our dynamical simula- 
tion to study the quenching effect is reported in 
refs. 



2. Numerical Simulation 

In Table. ^ simulation parameters used in this 
calculation and preliminary results of basic quan- 
tities are summarized. Our strategy of gauge 
generation with a rather fine lattice spacing and 
well-distributed topological charge is discussed in 
ref. 0. Since quark mass rrifa is introduced as 
a parameter of the boundary condition in the 
fifth dimension in domain- wall QCD, localiza- 
tion of chiral modes on both domain- walls in the 
fifth dimension tends to fail for a heavy quark 
mass. However, our small lattice spacing made 
the value of m c a acceptable as a domain-wall 
fermion: m c a ~ 0.45. We observed that, around 
this value, the shape of the wave function in the 
fifth dimension was qualitatively similar to that 
with much smaller m / . The small value of resid- 
ual quark mass m lcs S; 0.3 MeV demonstrates the 
good chiral symmetry. Therefore, it is expected 
that at least for those operators that do not mix 
with lower dimensional ones, finite L s causes neg- 
ligible chiral symmetry breaking effects. 




3. Kaon B-parameter Bk 

Bk on the lattice can be obtained as a ratio of 
the matrix elements of Qas=2 = [^7^(1 — 75)d] 2 
and axial current Aj. Our result is shown in Fig.^ 
as a function of raja. The fit function used is 
Bk = £o[l + Cm/aln(m/a)] + ^m/a with the 
coefficient C taken from analytic result [7 . In 
this chiral fit, x 2 /dof = 0.76 is obtained. 

Taking only statistical errors into account, our 
preliminary results are summarized as 

B^ tt (m f a = 0.0155) = 0.552(11), 



oRI-MOM 
B K 



MS,NDR 
K 



B 



(ji = 2 GeV) = 0.542(11), 
(p = 2 GeV) = 0.549(11), 



Bk = 0.764(15). 



Table 1 

Simulation parameters and preliminary results of 
basic quantities. 



size: 


24 3 x 48 


DBW2: 


= 1.22, 


#sweeps 


5,000 (B K ), 10,000 (K -> 7T, meson) 


DWF: 


M 5 = 1.65 , L s = 10 


m/a 


0.008 - 0.040, in step of 0.008 


m c a 


0.08, 0.12, 0.16, 0.20, 0.30, 0.40, 0.50 


#configs. 


77 (B K ), 23 (NPR), 42 (m rcs ), 




50 (K — > 7r, meson) 


a- 1 


2.86(9) GeV (m/a = 0, from p-meson) 




9.73(4) • 10" 5 (m/a = 0) 



The lattice value _B^ U is quoted at m/ = m s /2 
(filled symbol in Fig. ^) . By a non-perturbative 
renormalization procedure following 0, we ob- 
tained renormalization factor Z|P~ MOM (^ = 
2 GeV) = 0.9816(79). In this step, we also es- 
timated the degree of the mixing of the operators 
with wrong chirality. As expected by the small 
value of TO res , all elements of the amputated four- 
point function in the chirality basis were observed 
to be less than 0.2% of the (VV + AA, VV + AA) 
element which corresponds to Zb k - Though we 
shifted Bf MOM to Bjf [Ej and the renormal- 
ization group invariant (RGI) value Bk |11| with 
Nf = 0, results will change less than 1% even 
with N f = 3. Our result of B^(^ = 2 GeV) 
is consistent with CP-PACS (a" 1 ~ 2.9 GeV, 
24 3 x 60 x 16) and the previous work of RBC g|. 



4. K — > 7r matrix elements 

At the lowest order of chiral perturbation the- 
ory, K — > 7T7T matrix elements are proportional 
to K — » 7r matrix elements calculated on the lat- 
tice [3]. For i = 1 - 6, 9 and 10, these matrix 
elements are related as 
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where Q| /_0 ' 2 ^ is the contribution to the final 
state with / = 0, 2. Only for the case of AI = 
1/2, or J = 0, subtraction is needed to resolve 
mixing with a lower dimensional operator: 
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In Fig. |21 K — > 7r matrix elements of (left) 
and Q 2c = (su) l{u<£) l (right) are presented be- 
fore and after the subtraction of (J2J. In partic- 
ular, one finds that the slope of the subtracted 
matrix element of Q 2c is much smaller than that 
of Q 2 , which might mean a minor contribution 
of (tttt \Q 2c \ K) to AI = 1/2 rule. 



Fig. [3] shows 



K 



without the sub- 



K°)/(0\s l5 d\K°) 
rrida (right) whose slope 



traction (left) and 

as a function of m s a 
should be ae . Though these are an example with 
m c a = 0.40, their dependence on m c a is not vis- 
ible, so far. While the slopes of both quantities 
are determined within the error of 2% and 0.2% 
respectively, K — > tttt matrix element which is 
obtained as a combination of them has the er- 
ror more than 200% due to a severe subtrac- 
tion. Therefore, before we can quote result for 

K ^ , we need to significantly im- 
prove our statistics. 



For i = 7 and 8, K - 
elements are related as 



7T7T and K — > 7r matrix 
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We show ( 7r + 



+0(ml s ). (5) 

in Fig.Q]and observe its 

expected shape and intercept roughly consistent 
with previous works 02]. 
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